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The effects of mass imbalance in a three-component Hubbard model are studied by the dynamical 
mean-field theory combined with exact diagonalization. The model describes a fermion-fermion 
mixture of two different particle species with a mass imbalance. One species is two-component 
fermion particles, and the other is single-component ones. The local interaction between particle 
species is considered isotropically. It is found that the mass imbalance can drive the mixture from 
insulator to metal. The insulator-metal transition is a species-selective-like transition of lighter mass 
particles and occurs only at commensurate particle densities and moderate local interactions. For 
weak and strong local interactions the mass imbalance does not change the ground state of the 
mixture. 
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I. INTRODUCTION 


One of the fascinating problems in condensed matter 
physics is the metal-insulator transition (MIT). Espe¬ 
cially, the MIT driven by electron correlations has at¬ 
tracted a lot of attention. The electron correlations 
suspend the double occupancy of electrons, and the 
suspension causes the electron localization [l|. With 
the achievements of ultracold atom techniques, optical 
lattices can experimentally simulate various theoretical 
models of electron correlations @]. They have been pro¬ 
viding a novel stage for studying correlation effects in 
materials. In particular, the Mott insulator was ob¬ 
served in optical lattices of fermionic 40 K atoms with 
two hyperfine states and repulsive interaction between 
them [j, [ij. These optical lattices really simulate the 
two-component Hubbard model and they provide a con¬ 
nection between experimental observations and theoret¬ 
ical predictions. The optical lattices can also be es¬ 
tablished with different particle species which can be 
extended to have both large hyperfine multiplets and 
mass imbalances. For example, a mixture of single-spin 
state 40 K immersed in two-component fermionic atoms 
6 Li or a mixture of two-component state 171 Yb and six- 
component state 173 Yb have already been achieved [1, 6). 
Such achievements provide possible realizations of the 
MIT in multi-component correlation systems. Theoreti¬ 
cal studies already predicted a MIT in three-component 
Hubbard models [ 3 , Q. The MIT occurs at commen¬ 
surate particle densities when the local interaction is 
isotropic [7J. With an anisotropy of the local interac¬ 
tions, the MIT is also found at incommensurate half fill¬ 
ing [8J. However, in these studies, all component par¬ 
ticles have the same mass. Experiments can also sepa¬ 
rately tune the individual effective mass of each particle 
species and establish imbalanced mass mixtures. Indeed, 
in a mixture of 6 Li and 40 K atoms, the mass imbalance 
can be tuned in a wide range Q. The mixtures of two- 
component particles with different masses in the optical 


lattice generally lead to the difference between the hop¬ 
ping amplitudes associated with each component. With 
deviation from the balance limitation, some phase tran¬ 
sitions might happen and change the ground state. Low- 
temperature properties of the optical lattice in the influ¬ 
ence of the imbalance therefore are important and need 
to be considered. Indeed, mass imbalance in the two- 
component Hubbard model has been studied in detail 
in the one-dimensional and three dimensional cases that 
clarified the competing of some order states such as su- 
perfluid, ch arg e-, and spin-ordered states in the optical 
lattices 10-tl3|. Increasing the hopping imbalance or in¬ 
creasing difference between the bandwidth of each parti¬ 
cle species might lead the systems to a situation that one 
species is transformed from the metallic to the insulating 
state due to electron-electron correlations while the other 
still remains in its metallic or insulating state. The tran¬ 
sition is similar to the orbital-selective MIT in correlated 
multiorbital systems EM2. However, in ultracold mul¬ 
ticomponent mixtures the MIT may deal with odd num¬ 
bers of components, whereas in multiorbital systems it is 
impossible 0, [§]. Therefore, studying the effects of mas- 
simbalance on the MIT therefore is an interesting subject 
of the optical lattice systems. Actually, the effects of the 
mass imbalance on the MIT in fermion mixtures of two 
single-component particle species have been studied Ed- 
[l3j . The two species show distinct properties, which de¬ 
viate from the behaviors of balanced mass mixtures E3- 
In fermion-fermion mixtures of multicomponent particle 
species with extreme mass imbalance, where one particle 
species is extremely heavy and immobile, distinct MITs 
were also found [181 ] . The MIT can occur at both com¬ 
mensurate and incommensurate particle densities. De¬ 
spite that one species is always localized, the MIT can 
be classified as a collective, species-selective or inverse 
transition, depending on the local interaction and parti¬ 
cle densities [18| . 

In this paper we study the effects of the mass imbalance 
on the MIT in multicomponent fermion-fermion mix¬ 
tures. The mixtures consist of two-component and single- 
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component fermion particles. In contrast to the previous 
study [3, both particle species have finite masses and 
mobilities. Such fermion-fermion mixtures can be real¬ 
ized by 40 K with 6 Li atoms, or light 6 Li or 40 K with heavy 
fermionic isotopes of Sr or Yb. To model these mixtures 
with a mass imbalance we propose a three-component 
Hubbard model with different hopping amplitudes. Ac¬ 
tually, the hopping amplitudes can experimentally be 
tuned by the lattice potential and the recoil energy 0- 
In the balanced mass mixtures, the local interaction can 
drive the mixtures from metal to insulator states at com¬ 
mensurate particle densities [?]. In the extreme case, 
where one species is completely localized or nonhopping, 
the Hubbard model above reduces to a three-component 
Falicov-Kimball model and the MIT is found existing in 
both states of commensurate and incommensurate parti¬ 
cle densities 0]. Studying the MIT in between, i.e., with 
the hopping imbalance in the systems, therefore seems 
to be important to release the low-temperature quantum 
properties under the competition between the complex of 
electronic kinetic energy and the strong correlations. 

We study the three-component Hubbard model by em¬ 
ploying dynamical mean-field theory (DMFT) with exact 
diagonalization (ED). The DMFT has been used suc¬ 
cessfully to study the strongly correlated electron sys¬ 
tems 0, Hl| . The previous studies of MIT in the three- 
component Hubbard and Falicov-Kimball models were 
also based on the DMFT [3, B 0] ■ The DMFT is exact 
in infinite dimensions and fully captures local dynamical 
fluctuations 0,0]. However, it loses nonlocal correla¬ 
tions in finite dimensions. Within the DMFT, we found 
a MIT that is solely driven by the mass imbalance. This 
transition is a species-selective-like transition of lighter 
particles and occurs only at commensurate total particle 
densities and moderate local interactions. For weak and 
strong correlations, the mass imbalance cannot drive the 
mixture out of its state. 

The present paper is organized as follows. In Sec. II 
we describe the three-component Hubbard model with a 
mass imbalance. In this section we also present the ap¬ 
plication of the DMFT to the proposed Hubbard model. 
Numerical results for the detected MIT is presented in 
Sec. III. Finally, the conclusion is presented in Sec. IV. 


II. THREE-COMPONENT HUBBARD MODEL 
AND ITS DYNAMICAL MEAN FIELD THEORY 

We consider a three-component Hubbard model, the 
Hamiltonian of which reads 

H ^ ' taCictCja V ~ ^ , (1) 

where c\ a (cj a ) is the creation (annihilation) operator for 
the fermionic particle with hyperfine multiplet a at site 
i. a takes three different values, for instance, a = 1,2, 3. 
n-i a = c\ a Ci a is the number operator of the a-component 


fermionic particles at site i. t a is the hopping parameter 
of the a-component fermionic particles. U is the local 
interaction between the three component states of parti¬ 
cles. A common chemical potential /x is also introduced 
to control the total particle density n = ’Y^ ia {ni a )/N, 
where N is the number of lattice sites. The three- 
component Hubbard model can be realized by loading ul¬ 
tracold fermionic atoms with three hyperfine multiplets 
or fermion-fermion mixtures of different atomic species 
into optical lattices. However, in the Hamiltonian in Eq. 
©, the trapping potentials in the optical lattices are not 
considered. 

The mass imbalance solely depends on the hopping pa¬ 
rameters t a . In the three-component Hubbard model, the 
mass imbalance actually means the difference of the hop¬ 
ping parameters. In optical lattices, the particle hopping 
is established by the particle tunneling between nearest 
neighbor lattice potential wells. It can be tuned by the 
lattice potential amplitude V)( at and the recoil energy E ra 
of each component state of particles 0 

t a K, -^=E ra vl /4 exp(-2^u^), (2) 

V7T 


where v a = V^ at /E ra . The recoil energy E ra = k 2 /2m a , 
where k is the wave number of the laser forming the op¬ 
tical lattice, and m a is the mass of the a- component 
particles. The lattice potential amplitude V^ at can be 
different for different hyperfine states of particles. As 
a result, the hopping parameters can also be different 
even for the hyperfine states of the same particles with 
identical masses. In the following, we will consider the 
hopping imbalance t\ ^ ti = t%. This case can be in¬ 
terpreted as a fermion-fermion mixture of two different 
particle species. One species is particles with a single hy¬ 
perfine state (a = 1), while the other is particles with two 
hyperfine states (a = 2,3). Such mixture can be realized 
by loading fermion atoms 40 K and 6 Li, or of light atoms 
6 Li or 40 K with heavy fermion isotopes of Sr or Yb, into 
optical lattices. We parametrize the hopping amplitudes 
by 


t\ + t2 


At = 


t2 ~ t\ 
t\ + t2 


( 3 ) 

( 4 ) 


t is the average hopping amplitude of two particle species, 
and At describes the mass imbalance between them. It 
is clear that — 1 < At < 1. At = ±1 are the extreme 
mass imbalance, where one particle species is extremely 
heavy and localized 0 ]. Actually, the effective mass of 
the particles is inversely proportional to their hopping 
parameter. At > 0 indicates that the two-component 
particles are lighter than the single-component particles, 
and vice versa for At < 0. At = 0 is the balanced mass 
mixture and the model corresponds to the isotropic three- 
component Hubbard model [7]. In the balanced mass 
case, with sufficiently strong local interaction, the Mott 
insulator exists at commensurate particle densities [7j]. 
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The mass imbalance parameter At can experimentally 
be tuned in a wide range. For instance, in a mixture of 
6 Li and 40 K atoms, At can vary from 0.3 to 0.85 0. 

We study the three-component Hubbard model by em¬ 
ploying the DMFT. Within the DMFT the self-energy is 
a local function of frequency. It is exact in infinite dimen¬ 
sions. However, in finite dimensions, the DMFT neglects 
nonlocal correlations. The DMFT is well described in 
the literature, for example, in Ref. dj], For self-contained 
purposes, we present here a description of applying the 
DMFT to the three-component Hubbard model. The 
Green’s function of the a-component particles reads 


G a (k, iuj n ) 


1 


(5) 


where uj n is the Matsubara frequency, £k = 
S < ij > ex P(*k • (r, — i-j)) is the lattice structure fac¬ 
tor, and H a (iuj n ) is the self-energy. The local self¬ 
energy is determined from the dynamics of a single three- 
component particle embedded in a dynamical mean field. 
The action of the effective single particle reads 


P P 

Simp = - J J drdr' - T , )c a (r , ) 

oo a 

P 

jj r _ 

+— / d-r n a( T )n a '(T), (6) 

o a * a ' 

where Q a (r) is a Green’s function which represents the 
dynamical mean field. Q a (t) plays as the bare Green’s 
function in relation to the local Green function. It relates 
to the self-energy and the local Green’s function by the 
Dyson equation 


in Eq. ([6]) be identical to the local Green’s function in 
Eq. ©; i.e., 


G™ p (iuj n ) = G„(*w n ). (11) 


This equation completes the set of self-consistent equa¬ 
tions for the Green’s function. It can be solved numer¬ 
ically by iterations [2lfl . The most time-consuming part 
is the solving of the action <Si mp in Eq. ©• There are 


action S m 


21L 

culate it 

21 

, 22] 


The action in Eq. © is essentially 


equivalent to the Anderson impurity model [2ll . 


Haim = -p E c « c « + Tf E 


TlryTlcv' 


ct^ot' 


^ ^ ^pot^pcy^ot. H.C. + ^ ^ Ep 


p,a 


p,a 


( 12 ) 


where aj, a ( a pa ) is the creation (annihilation) operator 
which represents a conduction bath with energy level 
E pa - V pa is the coupling of the conduction bath with 
the impurity. The connection between the Anderson im¬ 
purity model in Eq. (fl2l) and the action in Eq. © is the 
following identity relation of the bath parameterspll. |22| 


E 



XUJn E pa 


Aa(A n ), 


(13) 


where \ a ( iu> n ) = ico n + p — Q ~ 1 {iui n ). The ED limits the 
conduction bath to finite n s — 1 orbits (p = 1, 2,..., n s — 1). 
Then A (iui n ) is approximated by 


G^iiuin) = G^iiuin) + T, a (iuj n ). (7) 

Here, the local Green’s function is 


n s — 1 

A %‘XiUn) = £ 

P= 1 


I Apq ] 2 

Epa 


(14) 


Get 


J d£p 0 (s) 


iu) n + p 


1 

T, a (iuj n ) + t a e' 


( 8 ) 


The bath parameters are determined from minimization 
of the distance d between A a (ioj n ) and A ^ s \iuj n ), 


where po(e) = J] k <5(£ — £ k) the bare density of states 
(DOS). Without loss of generality, we use the semicircu¬ 
lar DOS 


1 

d = WTl E < fc |A«(^n) - A^(* w „)| 2 , (15) 

' 1 n —0 


Po(e) = — \/l — £ 2 - (9) 

7r 

With the semicircular DOS, from Eqs. © and © we 
immediately obtain the Green’s function representing the 
dynamical mean field Q(iu} n ) from the local Green’s func¬ 
tion [2lJ 


Qa (Ay,,) — A n T p ^ G a (A n ). 


( 10 ) 


The self-consistent condition of the DMFT requires that 
the Green’s function obtained from the effective action 


where M is a large upper cutoff of the Matsubara fre¬ 
quencies dH HI- The parameter k is introduced to 
improve the minimization at low Matsubara frequen¬ 
cies. In particular, we take k = 1 in the numerical cal¬ 
culations. When the bath parameters are determined, 
we calculate the Green’s function of the Anderson im¬ 
purity model in Eq. (fl2l) by ED dH [H[. We also 
calculate the inter-species double occupancy Di nte r = 
( n n( n i 2 +n-i 3 ))/N, and the intra-species double oc¬ 
cupancy Z?; ntra = (n^n^/N. These double occupan¬ 

cies show the rate of the lattice sites occupied by two 
particles, and are experimentally accessible [3|. 
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FIG. 1. (Color online) The total particle density n, calculated 
by DMFT+ED (lines) and by DMFT+QMC (symbols), as a 
function of the chemical potential fi for different values of 
interaction U in the balanced mass case At = 0 (T = 0.025t). 
The DMFT+QMC results are reproduced from Ref. @. 


III. METAL-INSULATOR TRANSITION 

In this section we present numerical results analyzing 
the MIT under the influence of the mass-imbalance in¬ 
dicated by At [c.f. Eq. ©]■ The numerical results are 
obtained by the DMFT+ED with 3 bath orbits per one 
impurity component (i.e., n s = 4). Actually, the compu¬ 
tational time grows quickly with n s , since the impurity 
has 3 components. For the single-band Hubbard model, 
the DMFT+ED shows that two bath levels per one im¬ 
purity component usually give adequate results [23(. To 
check the accuracy of the ED, we compare our results 
with the ones obtained by the DMFT plus the Hirsch- 



FIG. 2. (Color online) The double occupancy as a function 
of the chemical potential /r, calculated by DMFT+ED (lines) 
and by DMFT+QMC (symbols) for different values of inter¬ 
action U in the balanced mass case At = 0 (T = 0.025t). The 
DMFT+QMC results are reproduced from Ref. j7. 


Fye quantum Monte Carlo (QMC) simulations in the 
case of mass balance (At = 0) [zj. The QMC simula¬ 
tions are an exact impurity solver, and unlike the ED 
they do not suffer the finite-size effects [2l|. First, Fig. [Q 
illustrates the comparison of the total particle density 
n as a function of the chemical potential /i obtained by 
DMFT combining with ED and QMC for two different 
U values. When the line n(/z) exhibits a plateau, it in¬ 
dicates an insulating state. The width of the plateau is 
equal to the insulating gap. Figure [T] apparently shows 
an excellent agreement between the ED and QMC simu¬ 
lations for both metallic and insulating phases in a whole 
range of the total particle density. In Fig. [2] we continue 
specifying the efficiency of the DMFT+ED calculation to 
consider the MIT in the Hubbard model by putting its 
beside DMFT+QMC results of the intra-species double 
occupancy Z?; n t ra . As a function of the chemical potential 
lb .Dintra again recovers the excellent agreement between 
the ED and QMC simulations for both metallic and insu¬ 
lating phases. The double occupancy is suppressed in the 
insulating phase. However, it still remains finite. The fi¬ 
nite value of the double occupancy in the insulating phase 
therefore is not a finite-size effect of the ED impurity 
solver. As pointed out in the literature the Brinkman- 
Rice approximation shows that_the double occupancy 
vanishes in the Mott insulator [24j. In reality a finite lo¬ 
cal interaction always allows virtual hoppings, that pro¬ 
duce very small but nonzero double occupancy in the 
insulating state. That feature has been addressed by the 
DMFT [Hi, I 25 U 27 I ] . The double occupancy vanishes only 
at the strong-interaction limit. Actually, the Brinkman- 
Rice approximation is based on the Gutzwiller variational 
wave function and the Gutzwiller approximation for eval- 



co/t 


FIG. 3. (Color online) The density of states (DOS), calcu¬ 
lated by DMFT+ED (lines) and by DMFT+QMC (symbols) 
for different chemical potentials fi in the balanced mass case 
At = 0 (17 = 3f, T = 0.05t). The DMFT+QMC results are 
reproduced from Ref. 0. 
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FIG. 4. (Color online) The total particle density n as a func¬ 
tion of the chemical potential /i for positive mass imbalances 
(At > 0) at fixed interaction U = 4t (T = 0.02f). 


uating the variational ground-state energy, and it admits 
the vanishing of the double occupancy in the insulating 
phase [2J, [29}. However, at finite dimensions, without 
the additional Gutzwiller approximation, the Gutzwiller 
variational wave function always produces a finite dou¬ 
ble occupancy for any local interaction [30{ . Moreover, 
extending DMFT for the finite-dimension case where the 
nonlocal correlations are taken into account has also il¬ 
lustrated the incomplete suppression of the double oc¬ 
cupancy in the Mott insulator (3lL i+1. In experiment, 
suppression of the double occupancy in the MIT of the ul¬ 
tracold two-component fermion atoms has been observed, 
but a small number of lattice sites, typically a few per¬ 
cent, still remains doubly occupied in the Mott insula¬ 
tor Q- To complete the comparison, Fig. [3] illustrates 
the density of states (DOS) of the particle component 
ImG 0 (« — ir))/n, evaluated by both DMFT+QMC Q 
and DMFT+ED. In the DMFT+ED calculation, we used 
77 = O.Olf. Despite the spiky structure in the DOS, the 
main features of the DOS obtained by the DMFT+ED 
resemble the ones shown by the DMFT+QMC. Inspect¬ 
ing the DOS at the chemical potential (i.e., the DOS at 
ui = 0) we see that the quantities determined by both ED 
and QMC simulations fit well together. In the metallic 
phase the DOS at the chemical potential level is always 
finite, while it vanishes in the insulating phase. From the 
above comparisons of ED and QMC results, we conclude 
that the ED solving for the effective impurity problem 
with n s = 4 in the DMFT already gives adequate results. 
In the following, we extend the calculation to the mass 
imbalance case (i.e. At 7 ^ 0) and then discuss its effects 
on the MIT in the three-component Hubbard model. 

In Fig. [1] we plot the total particle density as a function 
of the chemical potential when the mass imbalance varies 
in the region At > 0 at a given isotropic local interac¬ 
tion U. The value of the local interaction is chosen that 
allows a MIT. In the case of U = 4t we observe that the 


FIG. 5. (Color online) The total particle density n as a func¬ 
tion of the chemical potential /i for negative mass imbalances 
(At < 0) at fixed interaction U = 2.5f (T = 0.02t). 


mass imbalance possibly drives the mixture only from 
insulator to metal. Figure [4] also shows that when the 
mass imbalance is absent, At = 0, the line n(fi) exhibits 
plateaus at n = 1 and n = 2. The mixture is in the 
insulating phase at the commensurate densities only [ 7 }. 
With increasing the mass imbalance At, the plateaus re¬ 
duce and then disappear at large mass imbalances. These 
behaviors show a transition from insulator to metal solely 
driven by the mass imbalance at the commensurate den¬ 
sities. Actually, when At increases, the two-component 
particles become lighter and the single-component parti¬ 
cles become heavier. Due to this property, the mobility 
of the two-component particles tends to increase too. At 
large mass imbalances it overcomes the local interaction, 
and the mixture becomes metallic. In the At = 1 limit, 
the single-component particles are localized due to the 
vanishing of their hopping, but the two-component par¬ 
ticles can be in the metallic phase 0. In such way, the 
mass imbalance can only drive the mixture from insulator 
to metal. It cannot drive the mixture out of the metallic 
state. 

Next, we continue considering a dependence of the to¬ 
tal particle density on chemical potential but in the op¬ 
posite region of the mass imbalance (At < 0). Figured] 
shows us that for a given local interaction, U = 2.5 1, 
the mass imbalance can only drive the mixture from 
metal to insulator. With the increasing of |Af|, the 
two-component particles become heavier, and the single¬ 
component particles become lighter. When |At| remains 
being small, n(/r) does not show any plateau and the 
system is in the metallic state. The plateaus indicating 
an insulating state only appear if |At| is larger than a 
critical value. That happens at the commensurate den¬ 
sities n = 1 and n = 2, similar to the opposite situation 
with At > 0 (see Fig. |4}. In the limit of At = —1, 
the two-component particles become localized due to 
the vanishing of their hopping. However, the hopping 
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FIG. 6. (Color online) The density of states (DOS) of the two- 
component particles (blue solid lines) and of single-component 
particles (red dotted lines) for positive mass imbalances (At > 
0) at fixed total particle density n = 1 and interaction U = 4f 
(T = 0.021). 


of single-component particles remains finite. For weak 
interspecies interactions, the single-component particles 
are in a metallic phase. However, with increasing the in¬ 
terspecies interaction, the band of the single-component 
particles is split into two subbands, separated by a gap. 
This constitutes a MIT driven by inter-species particle 
correlations. Actually, the model with At = — 1 is equiv¬ 
alent to the spinless Falicov-Kimball model for strong 
interactions [33|. The spinless Falicov-Kimball model 
also exhibits a Mott-like MIT by splitting the conduc¬ 
tion band into two subbands [33 h35I |. 

The MIT can be understood in more detail in a feature 
of the DOS. In Figs. [G] and [7] we plot the DOS of parti¬ 
cles at the commensurate density n = 1. The other com¬ 
mensurate density n = 2 is equivalent to n = 1 via the 
particle-hole symmetry. Figures [6] and [7] confirm again 
the MIT driven by the mass imbalance in both regions 
At > 0 and At < 0. When the mass imbalance is ab¬ 
sent, At = 0, particles of two species equally participate 
in the MIT driven by the local interaction [7[. However, 
as addressed before in Fig. [I] and Fig. [5l in the region 
of At > 0, the mass imbalance can drive the mixture 
only from insulator to metal, whereas in the region of 
At < 0, it only drives the mixture from metal to insula¬ 
tor. In the At > 0 region, with increasing At the two- 
component particles become lighter, and the bands are 
mostly occupied by them. This indicates that the mass 



FIG. 7. (Color online) The density of states (DOS) of 
the two-component particles (blue solid lines) and of single¬ 
component particles (red dotted lines) for negative mass im¬ 
balances (At < 0) at fixed total particle density n = 1 and 
interaction U = 2.5 1 (T = 0.02t). 


imbalance also induces a population imbalance between 
the two particle species. In the extreme limit, At = 1, 
the single-component particles are localized, and only the 
two-component particles take part in the MIT jig}. We 
can refer to the MIT driven by the mass imbalance as 
a species-selective-like MIT. In this MIT the lighter par¬ 
ticles are dominantly involved in driving the mixtures 
to the metallic state. However, this species-selective-like 
MIT is different in comparison with the orbital-selective 
MIT [HjT7t . In the orbital-selective MIT, the wide band 
remains metallic when the narrow band becomes insulat¬ 
ing. In this selective-species-like MIT, the bands of both 
species are insulating and when the MIT occurs, they 
together become metallic, but the lighter particles are 
dominant. Only at the limit At = 1, the heavier parti¬ 
cles are always localized. 

Similarly to the region of At > 0, when At < 0, the 
mass imbalance also induces a population imbalance be¬ 
tween two particle species, but this population imbal¬ 
ance is not strong as in the region of At > 0. The band 
filling of the single-component particles, which become 
lighter with increasing |At|, is larger than the one of the 
two-component particles. At the limit of At = —1, the 
two-component particles are localized due to vanishing 
of their hopping amplitude. We also refer to the MIT 
driven by the mass imbalance in the At < 0 region as 
a species-selective-like MIT, because the lighter particles 
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are dominant over the heavier ones in the transition. In 
contrast to the orbital-selective MIT, where particles of 
both spins always participate in the transition HM3, 
in this At < 0 region, particles of a single component 
are dominantly involved in the MIT. It is similar to the 
Mott-like MIT of the spinless electrons in the Falicov- 
Kimball model [33j - l35j . The MIT occurs when the band 
of the single-component particles is split by a gap due to 
the interspecies interaction. 

Experiments could observe the mass-imbalance-driven 
MIT by measuring the doubly occupied sites. Indeed, the 
MIT of ultracold 40 K atoms was detected by measuring 
the double occupancy Q. In Fig. [8] we plot the intra- and 
interspecies double occupancies at total particle density 
n = 1 in both At > 0 and At < 0 regions. One can 
see a kink of these double occupancies at the point of 
the MIT. It indicates that the MIT is a first-order phase 
transition. In the region of At > 0, the right panel of 
Fig. [8] shows us that the double occupancies in the in¬ 
sulating phase are small, but nonzero, as we have dis¬ 
cussed previously. With increasing the mass imbalance, 
at the MIT, the intraspecies double occupancy abruptly 
increases, while the interspecies double occupancy de¬ 
creases to zero. Actually, with large mass imbalances 
the metallic state is dominantly occupied by the two- 
component particles. As a consequence, the intraspecies 
double occupancy increases, while the interspecies dou¬ 
ble occupancy tends to vanish when At —> 1. In contrast 
to the region of At > 0, in the region of At < 0 (the 
left panel of Fig. |SJ), the intraspecies double occupancy 
decreases to zero, when | At| increases. It is clear that in 
this region, the two-component particles become heavier, 
and in the At —> — 1 limit, they are actually localized. 



At 


FIG. 8. (Color online)The intraspecies double occupancy 
Sintra (red filled circles) and the interspecies double occu¬ 
pancy Sinter (blue filled squares) at fixed total density n = 1 
(T = 0.02t). The left panel plots the region of At < 0 
(U — 2.51), while the right panel plots the region of At > 0 
(U = 4 1). 


Consequently, the intraspecies double occupancy tends 
to vanish when At —» —1. However, the interspecies 
double occupancy remains finite in the insulating phase. 
Actually, in the region of At < 0, the mass imbalance 
drives the mixture from the metallic to the insulating 
state. Therefore, the MIT may occur only for the inter¬ 
mediate local interaction that should be smaller than the 
critical value of the local interaction for the MIT in the 
mass balanced case. The local interaction in such value 
range is not strong enough to suppress the double occu¬ 
pation. Instead, it allows virtual hoppings that produce 
double occupation of interspecies particles, despite that 
the band of the single-component particles already opens 
a gap. 

To characterize the influence of the correlation-driven 
MIT in the case of mass imbalance, in Fig. [9] we show 
the dependence of the double occupancies on the lo¬ 
cal interaction for different At. In the balanced mass 
case (At = 0), Tenter = 2Di ntra - Both inter- and intra¬ 
species double occupancies abruptly change at the MIT, 
and they are suppressed in the insulating state. In the 
presence of the mass imbalance, only the intra- or inter¬ 
species double occupancy exhibits a kink at the MIT. In 
the region of At > 0, the double occupancy of the two- 
component particles also indicates the rapid suppression 
at the MIT, while the interspecies double occupancy con¬ 
tinuously decreases with increasing the local interaction. 



FIG. 9. (Color online) The intraspecies double occupancy 
73intra (red filled circles) and the interspecies double occu¬ 
pancy -Dinter (blue filled squares) at fixed total density n = 1 
and different mass imbalances (T = 0.02t). 
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FIG. 10. Phase diagram at total particle density n = 1 (T = 
0 . 021 ). 


These behaviors interchange each other in the opposite 
region with At < 0. They indicate the active role of 
the two-component particles in the At > 0 region and of 
the single-component particles in the region of At < 0 
in the MIT. It means that only lighter particles are ac¬ 
tively involved in the MIT. In the insulating phase, the 
double occupancies are suppressed, but they remain fi¬ 
nite and only vanish at the strong interaction limit. The 
double occupancies are experimentally accessible, but it 
is a challenge to experimental observations of the kink of 
the double occupancies. 

In Fig. [10] we summarize a phase diagram expressing 
the MIT in the (At, U ) plane at the commensurate to¬ 
tal particle densities. The phase diagram is constructed 
from the dependence of the total particle density on the 
chemical potential as discussed before. The insulating 
phase is detected when a plateau appears in the line of 
n(/Lt). In contrast to the two-component Hubbard model, 
we do not observe any phase coexistence of metal and 
insulator [26], [28[ . By increasing At the critical U value 
monotonically increases. However, since — 1 < At < 1, 
the MIT driven by the mass imbalance can occur only at 
a finite moderate range of the local interaction. For weak 
and strong local interactions, the mass imbalance cannot 
drive the mixture out of its ground state. In fermion 
mixtures of 40 K and 6 Li atoms, At can vary from 0.3 to 
0.85 @]. One may expect to detect the mass-imbalance 


driven MIT in such ultracold atom mixtures at commen¬ 
surate particle densities and moderate local interactions. 


IV. CONCLUSION 

We have studied the MIT driven by the mass imbal¬ 
ance in the three-component Hubbard model. Within 
the DMFT with exact diagonalization, we have found the 
MIT driven by the mass imbalance at commensurate to¬ 
tal densities, like in the balanced three-component Hub¬ 
bard model |7|. The MIT is solely driven by the mass im¬ 
balance. The positive mass imbalance can only drive the 
mixture from insulator to metal, while the negative one 
drives the mixture from metal to insulator. In order to 
explore the MIT in the three-component Hubbard model 
with mass imbalance, we have also calculated the dou¬ 
ble occupancies of both two-component particles (Di ntra ) 
and different species particles (fainter) as functions of At 
and local interaction U. At a critical U value, the dou¬ 
ble occupancies exhibit a kink, indicating the MIT tran¬ 
sition of the first order. The more enhancing the mass 
imbalance, the MIT takes place at larger critical U value. 
The phase diagram expressing the MIT in the (At, U ) 
plane is also constructed. It shows that the MIT occurs 
only at moderate local interactions. For weak and strong 
local interactions, the mass imbalance cannot drive the 
system out of its ground state. Actually, the mass im¬ 
balance induces a population imbalance between particle 
species, and the lighter particles dominantly take part in 
the MIT. The MIT can be interpreted as a light-particle 
species selective transition. These features are distinct 
from the behaviors of balanced systems. We predict the 
MIT driven by the mass imbalance in fermion-fermion 
mixtures of ultracold atoms, loading into optical latices, 
for instance, the mixture of 40 K and 6 Li atoms. By mea¬ 
suring the ratio of the doubly occupied sites in the optical 
lattices, the MIT could be observed by tuning the mass 
imbalance ratio. 
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